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In early 2014, construction of a new watt balance, named NIST-4, has started at the National 
Institute of Standards and Technology (NIST). In a watt balance, the gravitational force of an 
unknown mass is compensated by an electromagnetic force produced by a coil in a magnet system. 

The electromagnetic force depends on the current in the coil and the magnetic flux integral. Most 
watt balances feature an additional calibration mode, referred to as velocity mode, which allows one 
to measure the magnetic flux integral to high precision. In this article we describe first measurements 
of the flux integral in the new watt balance. We introduce measurement and data analysis techniques 
to assess the quality of the measurements and the adverse effects of vibrations on the instrument. 


I. INTRODUCTION 


Watt balanceJIl have a long history at the National In¬ 
stitute of Standards and Technology (NIST). In 2011, it 
was decided to build a new watt balance (NIST-4)P1 which 
will be used to realize the kilogram once the international 
community has agreed to proceed with the redefinition of 
the international system of units, the Systeme Interna¬ 
tional d’Unites (SI)^. In one aspect, the new watt bal¬ 
ance is a departure from previous instrument^^H^ built at 
NIST because a permanent magnet system instead of an 
electromagnet is used to generate the required magnetic 
field. This design choice is also used in other watt bal- 
ance^^K^ and leads to an instrument that is much smaller 
than its predecessors. 

In the watt balance experiment, the gravitational force, 
F = mg, on a mass, m, at the local gravitational acceler¬ 
ation, g, is counteracted by an electromagnetic force pro¬ 
duced by a current carrying coil immersed in a magnetic 
field. The force generated by this coil can be expressed, 
in its simplest form, as F = IIB, where / denotes the cur¬ 
rent in the coil, B the average value of the magnetic flux 
density perpendicular to the wire and I the wire length 
of the coil. 

Balancing the gravitational with an electromagnetic 
force is one of two measurements, sometimes called 
modes, that are performed with a watt balance. The 
purpose of the second mode is to precisely measure the 
flux integral IB. This measurement is implemented by 
moving the coil perpendicular through the magnetic field 
while recording the coil’s velocity, v, and the induced 
voltage, V, in the coil. The ratio gives the flux integral, 
IB = Vjv. 

Construction of NIST-4 started after three years of de¬ 
sign and preparation in January of 2014. Immediately af¬ 
ter the balance was assembled, measurements of the flux 
integral were taken with the goal to demonstrate that 
a measurement of the flux integral with a relative un¬ 
certainty of 10“^ is possible. The main concern for this 
measurement was the effect of mechanical vibrations of 
the instrument on the data. NIST-4 is smaller and lighter 
than NIST-SP leading to higher frequencies of its mechan¬ 
ical resonances. Preliminary measurement of the ground 
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FIG. I. Three dimensional drawing of the NIST-4 watt bal¬ 
ance. For clarity the vacuum vessel and elements of the sup¬ 
port structure have been omitted. 


vibration at the NIST-4 site indicated larger amplitudes 
than at the NIST-3 site. In particular, the amplitude of 
the vertical ground motion at 29 Hz at the NIST-4 site is 
about 30 times larger than at the NIST-3 site, see Fig.[^ 
This vibration is most likely caused by electrical motors 
in large air conditioning units nearby. 


II. THE APPARATUS 


Figure shows a drawing of NIST-4. Identical to the 
previous three instruments, a wheel serves as the balanc- 
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ing mechanism. In the current setup, an aluminum wheel 
with a diameter of 60.96 cm and a thickness of 2.54 cm 
is used. A polished stainless steel ring is press-fitted on 
the outside diameter of the aluminum wheel to provide 
a hard, polished surface for a multi-filament band to roll 
on. The multi-filament band is made from 65 Platinum- 
Tungsten wires. The wheel is supported at its center 
from a 7.87 cm long tungsten-carbide knife edge that can 
pivot on a polished tungsten-carbide block. The center 
of mass of the wheel is adjusted to be closely beneath the 
knife edge. The free period of oscillation about the knife 
edge is about 10 s. 

The two sides of the wheel are referred to as the main- 
mass and the counter-mass side. On the main-mass side, 
the multi-filament band connects to a rotational coupler 
which is a bundle of 150 straight Platinum-Tungsten 
wires each with a diameter of 76 fim and an overall length 
of 89 mm. The coupler terminates at a flexure system 
comprised of two independent flexures that pivot about 
the same point. One flexure supports the mass pan, while 
the other flexure supports the spider, a three-pointed star 
made from aluminum. The spider connects to the coil 
by three carbon fiber rods that are mounted with two 
flexures on either end allowing the coil to move parallel 
to the spider. 

When the wheel is rotating, the multi-filament band 
rolls on or off the stainless steel surface shortening or 
lengthening the free hanging part of the band. Even 
though the band was annealed, internal stresses produce 
torques about the vertical axis on the lower end of the 
band. If not counteracted, these torques would produce 
an unwanted rotational oscillation of the coil about the 
z-diX.is. An electrostatic torque compensation system is 
used to mitigate this problem. The actuator consists 
of six high voltage electrodes, three for clockwise and 
three for counter clockwise torques, surrounding three 
grounded gold plated glass plates mounted on the carbon 
fiber rods. The input of the PID loop of the torque com¬ 
pensation system is calculated from two optical sensors, 
each monitoring the x and y position of a retro-reflector 
mounted on the coil former. From these data, the rota¬ 
tion of the coil around the ^ axis can be calculated. Such 
a system was already used successfully in the NIST-3 ex¬ 
periment. When the PID loop is off, the amplitude of the 
rotation can be adjusted with the rotational coupler to 
within 0.1 mrad when the coil moves up and down. With 
the PID loop on, the amplitude is better than 0.01 mrad. 
The rotational coupler is an element with low torsional 
stiffness that allows us to servo the rotational orientation 
of the coil with moderate voltages (below 1000 V). 

The coil former is made from glass-epoxy laminate, 
grade FR-4, the material from which printed circuit 
boards are typically made. The coil former contains two 
rectangular grooves, vertically displaced by 11.5 mm. In 
each groove, a coil with 464 turns and a mean radius of 
21.75 cm is wound, using AWG-24 wire (diameter approx¬ 
imately 0.5 mm), yielding a resistance of about 56.2^4. 
Electrically both coils are connected in series. 



FIG. 2. Topology of the heterodyne interferometer used to 
measure the vertical position and velocity of the coil. The 
dotted box denotes the vacuum envelope. Legend: BS = 
beam-splitter, PBS = polarizing beam splitter, HWP = half¬ 
wave plate, QWP = quarter wave plate, CCR = corner cube 
retro-reflector. 


The coil is inside a permanent magnet system, where 
two Samarium-Cobalt, Sm 2 Coi 7 , disks source the mag¬ 
netic field and soft iron (AISI 1021) guides the flux 
through a 3 cm wide annular gap. The magnetic flux 
density in the center of the 15 cm high gap is 0.55 T. 
The mass of the magnet system is approximately 816 kg. 
More details on the design and construction of the perma¬ 
nent magnet system can be found in referenceJI^ ancP^. 

From the coil, three aluminum posts extend down to 
a second spider which supports a corner cube reflector 
in its center. This reflector is used to interferometrically 
measure the coil’s vertical position and velocity. The 
light beam of the interferometer is traveling through a 
central bore in the magnet. 

On the counter-mass side of the wheel, a smaller coil 
is suspended in a magnet system to provide a small drive 
force for the velocity control. This coil has a nominal di¬ 
ameter of 18.3 cm 297 turns, with a resistance of 21 D. 
The coil is inside a smaller magnet system that uses 
Neodymium-Iron-Boron, Nd 2 Fei 4 B, permanent magnets 
to produce a flux of 0.14 T in the center of a 3 cm wide 
gap. In the middle of the 10.16 cm high gap, the flux 
integral is 23.6 T m. This counter-mass coil is connected 
by three rods to an aluminum cylinder that hangs on a 
multi-filament band off the wheel. The aluminum cylin¬ 
der is necessary to counter-balance the heavier main mass 
assembly. 

The vertical position of the coil is measured with a het- 
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FIG. 3. Timing diagram for the time interval analyzer (TIA) 
and the three digital voltmeters used in the experiment. 

erodyne interferometer using an Agilent 5517C laseras 
the light source. The optical layout, shown in Fig ure 
was designed to minimize polarization mixing^^. The 
half-wave plate in the optical path before the polarizing 
beam splitter allows us to align the polarization direc¬ 
tions of the light beam to the polarization axes of the 
polarizing beam splitter cube. The additional polarizing 
beam splitter in the reference arm reduces the amount of 
horizontally polarized light leaking in the reference arm of 
the interferometer, reducing the error produced by polar¬ 
ization mixing. The effect of this additional beam splitter 
can be clearly observed in the amplitude spectral density 
of the measured velocity and reduces the amplitude spec¬ 
tral density by about a factor of two. 

The signals from two optical pickups, one for the ref¬ 
erence and one for the interferometer signal, are con¬ 
verted to an electrical signal with a photo receiver (Ag¬ 
ilent E1709A). The electrical signals are fed into a time 
interval analyzer (TIA) made by Brilliant Instruments, 
Model BI221. The time interval analyzer is triggered by 
the electrical pulses received from the reference. The sig¬ 
nal to arm the trigger is generated by a pulse generator 
using a frequency up to 300 Hz, see Fig. The trig¬ 
gers for three digital voltmeters (DVMs) that are used to 
sample the induced voltage in the coil are generated at 
the same frequency, but a time delay is adjusted to ob¬ 
tain simultaneous measurements of voltage and velocity. 
This delay was adjusted to within 500 ns by moving the 
coil in a sinusoidal waveform and analyzing the Lissajous 
figures obtained by plotting voltage versus velocity. 

To understand the vibrational motion of the coil with 
respect to the magnet system, a closer look at the suspen¬ 
sion of the watt balance is necessary. A cross-sectional 
side view and top view of the main suspension elements 
can be found in Figure]^ The main watt balance rests 


on two sand filled steel tubes (inner diameter 17.8 cm, 
outer diameter 20.3 cm) which are supported by four legs. 
These legs are located on the corners of a rectangle with 
a length of 2.33 m and width of 0.81 m. The magnet en¬ 
gages in a pseudo-kinematic mount onto the beams. Two 
support points are mounted at the top diameter of the 
magnet, the third on a panhandle that rests on a bridge 
between the two beams. A large aluminum plate in the 
shape of a teardrop is bolted to the top surface of the 
magnet. This aluminum plate is the base for mounting 
the optics and is the building platform for the mass lift 
(not shown in Fig. . Three posts, again sand filled, 
stand on this plate and support the table which is a 
trapezoidal shaped plate with a central slot to accommo¬ 
date the wheel. On top of the table, rests the extremely 
large flexure (ELF). The ELF allows us to move the knife 
edge support by ±0.5 mm along the x- and y-direction. 
These degrees of freedom are necessary to align the coil 
concentric with the magnetic field in order to reduce the 
torques on the coil during the weighing mode. A detailed 
description of the ELF can be found irP^. 

The main vibrational frequencies in the balance listed 
in order of decreasing frequencies are: 

• 30 Hz. The suspension of the coil on the main mass 
side, consisting of the bands, the spider, and the 
carbon fiber rods has a bounce frequency of 30 Hz. 

• 21 Hz and its harmonics. The moving part (drawn 
in dark color in Fig. of the ELF is bouncing 
at this frequency. To damp this vibration the inner 
part has been pretensioned and damped with bolts. 
While this damps the 21 Hz vibration, it introduces 
vibrations at harmonics, particularly at 63 Hz. In 
the future, we plan to insert damping material be¬ 
tween the table and the ELF. 

• 16 Hz. The watt balance on the two steel beams 
moves up and down at a frequency of 16 Hz. 

• 1.7 Hz. This frequency is the double pendulum mo¬ 
tion from the point where the band rolls off from 
the wheel. In this mode, the spider moves in one 
direction and the coil in the other. 

• 0.85 Hz and 0.834 Hz. The tilt of the coil oscillates 
about its center of mass at two frequencies around 
the X- and y-axis, respectively. These motions are 
referred to as wobble motions of the coil. 

• 0.56 Hz. The pendulum motion of the coil, more 
precisely the shearing motion between the coil and 
the upper spider. For this motion, all six flexures 
on the carbon fiber rods are bending. 

The frequencies listed above were measured with the 
wheel at the middle position in its travel range. Rotating 
the wheel will slightly change some of these frequencies 
due to shortening or lengthening of the suspension. 
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FIG. 4. The support system of the watt balance. In the middle, top: A cross-sectional drawing as seen from the side. In the 
middle, bottom: A top view of the instrument. For simplicity the vacuum vessel is omitted. Other than through the four 
bellows shown, there is no connection between the balance and the vacuum vessel. The directions of our coordinate system are 
indicated in the graph. On the left in the dotted box: Enlarged truncated drawing of one of the three coil support rods. On 
the left in the dashed box: The termination of the multi-hlament band and rotational coupler. On the right: A top view of 
the extremely large flexure (ELF). The outer ring is bolted to the trapezoidal plate. The inner part of the ELF, drawn with 
darker color, can move in x- and y- direction by ±0.5 mm. 


III. DATA ANALYSIS AND FIRST RESULTS 

The induced voltage, V, in the coil is proportional to 
its vertical velocity, v, and the flux integral, IB. Because 
of the finite size of the magnet, the magnetic flux varies 
along the gap and thus the flux integral is a function 
of coil position. Hence, V{t^z) = v{t)lB{z). Neglect¬ 
ing the non-uniformity of the flux integral, during the 
watt experiment, the coil moves with a nominal velocity 
'^nom and a nominal induced voltage I4om ~ is 

measured. Here, IB is 697 T m and thus for a nominal 
induced voltage of 1V a velocity of 1.43 mm/s is required. 

For a stationary balance, i.e., a non-rotating wheel, the 
residual coil movement is a small vibration around an av¬ 
erage position 2 : 0 . In this case, the excursions from the 
nominal positions are small and the position dependence 
of IB can be neglected. However, the quotient of volt¬ 
age to velocity cannot be calculated because the velocity 
will sometimes be zero driving the quotient to infinity. 
Instead, we introduce a quantity, named the residual, 

e{t) = V{t)-v{t)lB{zo). (1) 


For a vibrationally driven coil and sufficiently long obser¬ 
vation time, the expectation value of e{t) is (^(t)) = 0. 
The lB{zo) can be calculated from the slope of the mea¬ 
sured voltage as a function of velocity, V(v). The residual 
is useful in two ways. First, in the frequency domain, it 
allows one to identify the regions where the motions of 
the coil track the induced voltage and where they don’t. 
Second, as we will show below, from measurements of 5 
an estimate of the uncertainty of the IB measurement in 
the real experiment can be calculated as follows. 

In velocity mode of a watt balance experiment, one 
tries to estimate the flux integral by measuring the in¬ 
duced voltage and the velocity of the coil, i.e. 


V 

IB=-. (2) 

V 

The relative uncertainty squared of the determination of 
the flux integral in the watt balance experiment can be 
written as 




( 3 ) 
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where denotes the correlation coefficient between the 
voltage measurement and the velocity measurement. The 
squared value of the absolute uncertainty of e is given by 


2 

al = (Jy + {IBY^-2lBryvcrvo-v (4) 

Dividing this equation by the squared value of nominally 
induced voltage during the watt experiment, 
using the relationship, Idiom = yields 
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The relative uncertainty achieved for measuring IB is the 
same as the quotient of the absolute uncertainty of 5 
divided by the nominal induced voltage. 



FIG. 5. Top Panel: Amplitude spectral density of the ground 
vibrations measured at the site of the new watt balance 
(NIST-4) and the old watt balance (NIST-3). Bottom Panel: 
Amplitude spectral density of the induced voltage V, the ve¬ 
locity multiplied with the flux integral, vlB, and the residual 
s. All three spectra are divided by Tnom = 1V. 


In order to investigate the basic resonances of the watt 
balance, a sinusoidal voltage with an amplitude of 0.5 V 
and a frequency sweeping from 1 Hz to 150 Hz was ap¬ 
plied to the counter mass coil, while the induced voltage 
in the main coil and the velocity of the main coil was 
recorded. Fig. shows the amplitude spectral densities 
of the voltage, the velocity, and the residuals with the 
same normalization as above. Two strong resonances 
can be seen, one at 30 Hz, the other at 67 Hz. Both 
resonances are double peaks. The 30 Hz is the up down 
vibrational mode of the coil caused by stretching of the 
support. The voltage and velocity are in phase and co¬ 
herent up to 80 Hz. At 80 Hz, the coherenc^I^ decreases, 
see upper graph in Fig. 



FIG. 6. Lower graph: The amplitude spectral density of the 
induced voltage, the coil’s velocity and the residual while the 
counter mass coil is driven by a swept sine wave with a fre¬ 
quency ranging from 1 Hz to 150 Hz and a 0.5 V amplitude. 
The normalization of the spectra is the same as in Fig. [^Up¬ 
per graph: coherence between voltage and velocity during this 
measurement. 


The amplitude spectral density of the three quantities, 
V, vlB^ and 5 are shown in figure All three spectral 
amplitudes are divided by Viom=l V. These three spec¬ 
tra correspond to the amplitude spectral density of the 
square roots of the term on left and the first two terms 
on the right of the equal sign in Eq. [4| All three spectra 
share the same vertical unit, 1 /VHz The traces for the 
spectra of V/Viom and vlB/Vnom are on top of each other 
for a large range in frequency. At these frequencies, both 
signals are in common mode and the amplitude spec¬ 
tral density of e/Vnom is much smaller than either one. 
With increasing frequency at about 8 Hz, the noise in 
the velocity measurement exceeds the noise in the volt¬ 
age measurement and dominates the amplitude spectral 
density, i.e., the amplitude spectral density of the resid¬ 
ual follows that of the velocity. In general, the amplitude 
spectral density of the velocity increases with increasing 
frequency, because the velocity is the derivative (multiply 
by / in frequency domain) of the position measurement, 
which exhibits a nearly flat amplitude spectral density. 


When the balance is servoed to a fixed position or 
moved with constant velocity, the twin peaks at 30 Hz 
dominate the amplitude spectral density of the residual 
and hence the amplitude spectral density of the flux in¬ 
tegral. These signals make it difficult to obtain a precise 
measurement of IB. Using longer integration times on 
the voltmeter and correspondingly larger sample times 
for the time interval analyzer did not solve the problem. 
In this case, the peaks are aliased into the measurement 
band introducing spurious signals. A promising solution 
is to sample the DVMs and the TIAs with a relatively 
large frequency (120 Hz) and to apply a digital, finite 
impulse response (FIR) low-pass filter to both channels 
before calculating the residuals or the ratio. Currently, 
we use a filter with a passband-edge frequency of about 
6 Hz and an attenuation of at least 110 dB for frequencies 
larger than a stopband-edge frequency of 10 Hz. 

To measure the flux integral, the coil is moved through 
the magnet. The balance is in servo control using the 
velocity mode motor as actuator. The set-point of the 
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position servo varies with a triangular waveform with an 
amplitude of 3 cm and a period of 96 s. To avoid bias 
at the turning points, the data analysis is limited to the 
inner 5 cm of travel, where the coil moves with a constant 
velocity of 1.36 mm/s. Figure shows the data of one 
such sweep. The IB is calculated by dividing the induced 
voltage by the measured velocity. In the top panel, the 
amplitude spectral density of the IB is shown for the 
filtered and unfiltered data stream. The amplitude in 
the unfiltered data is dominated by the amplitude density 
beneath the twin peaks at 29 Hz and elevated white noise 
from 10 Hz to the Nyquist frequency of 60 Hz. 

After filtering the data, three dominant peaks remain 
at 0.55 Hz, 0.85Hz and 1.7Hz, corresponding to the pen¬ 
dulum, wobble, and double-pendulum motions. These 
motions do not cancel in the ratio due to misalignment 
of the electrical center with the optical center. 

A common alignment practice is a two step procedure. 
In one step, the optical center is aligned to the mass cen¬ 
ter by observing coil wobble and minimizing the second 
harmonic in the optical signal. The difference between 
optical center and mass center is referred to as Abbe 
offset In another step, the electrical center is aligned to 
the mass center by injecting current in the coil and ob¬ 
serving the coil tilt. After the two step procedure, the 
electrical center is aligned with the optical center. In 
this preliminary setup, only a coarse alignment of the 
optical center with the mass center was possible. In the 
future, we plan to use three interferometers to track three 
points on the coil former, each pair spaced 120^ apart. 
By mathematically combining the readout of the three 
interferometers, the optical center is obtained. Changing 
the relative weighting of the three interferometer signals 
allows us to move the optical center in software, as re¬ 
quired by one of the steps above. 

In addition to reducing the Abbe offset, two other 
strategies to reduce the effect of parasitic motions on the 
measurement result have been used in the past: In NIST- 
3, active damping is used before and after each sweepl^. 
In the NRC balance, the starting time of the sweep is ran¬ 
domized with respect to the phase of the tilt motiorP^. 

From the filtered data, the flux integral as a function 
of coil position is calculated and shown in the lower plot 
of Figure The measured variation of the flux integral 
as a function of position agrees well with measurements 
performed while characterizing the magnet as described 
in RefP. The relative amplitude of the variation is 2 x 
10 -^. 

At low frequency, the amplitude spectral density of the 
flux integral is about I0“^Tm/\/Hz. Integrating this 
level over a 48 s long sweep yields a relative uncertainty 
in IB of 2 X 10 To verify this level of uncertainty, the 
coil was moved through the magnet for an entire night. 
A sixth degree polynomial lB{z) is adjusted in a least 
squares fashion to the data obtained in each sweep, see 
the black line in the lower plot of Fig.[^ This polynomial 
is then used to calculate IB{0). A plot of these numbers 
is shown in Figure The results for the coil moving 




z (cm) 


FIG. 7. Top graph: The amplitude spectral density of the 
measured IB during a constant velocity motion of the coil. 
The IB is calculated from the ratio of the measured voltage 
and velocity. In one case, dashed line, the ratio is calculated 
from the as-measured data stream. In the second case, solid 
line, the voltage and velocity data is passed through a digital 
FIR filter. The 3 dB corner frequency of the filter is at about 
6 Hz. Bottom graph: The flux integral calculated from the 
filtered data from one 5 cm long motion through the magnetic 
field. 

in one direction differ from that by moving in the other 
direction due to thermal voltages in the system. Ref.l^ 
shows the influence of thermal voltages in velocity mode. 
These thermal voltages, however, seem to remain fairly 
constant over the course of the night. Overall the flux in¬ 
tegral drifts by a few parts in 10“^. This drift is caused 
by temperature changes of the Sm 2 Coi 7 in the perma¬ 
nent magnet. These data were taken several days after 
the vacuum system had been pumped down, when the 
magnet has reached a fairly stable temperature. The rel¬ 
ative scatter in the data is close to the prediction from 
the amplitude spectral density, about 1.3 x 10“^. Hence, 
a measurement of IB with a relative statistical uncer¬ 
tainty of 1.3 X 10“^ should be possible in about 4,800 s 
at a velocity of 1.36 mm/s, assuming that the noise is 
white and stationary. At present, a relative statistical 
uncertainty of I x 10“^ would be achievable within 2.25 
hours of measurement time. 


IV. CONCLUSION 

The main components of the new NIST watt balance 
have been assembled. First measurements of the flux in¬ 
tegral were performed and the noise levels and structural 
resonances of the system were thoroughly investigated. 

We found that the coil motion has one strong vibra¬ 
tional resonance at 29 Hz. In this mode the coil moves up 
and down, stretching the bands and carbon fiber rods. 
Upon close inspection, both resonance peaks split into 
two peaks. These resonances dominate the ratio of volt 
to velocity, i.e., the flux integral. One successful way to 
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FIG. 8. Measurement of the flux integral throughout one 
night. 


remove excess noise from the data stream is to measure 
volt and velocity with a sample rate of 120 Hz and then 
use a digital low pass filter with a cut-off frequency of 
GHz. From the two filtered data channels, the ratio is 
then calculated. Since both channels are filtered through 
the same digital filter, no bias is added to the data. This 
procedure will, however, add a phase delay to the final 
values of 

In addition the performance of the system can be en¬ 
hanced by improving the alignment. Especially by using 
three interferometers to allow a alignment of the Abbe 
offset. 

With the current system, we have shown that it is pos¬ 
sible to measure the flux integral with a relative statisti¬ 
cal uncertainty of 1.3 x 10“^ for one sweep. 

In the future, we will implement the force mode and 
incorporate a Josephson voltage standard in the system. 
We hope to have a fully working watt balance by summer 
2015. 
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